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ABSTRACT 

We consider a class of spatially flat cold dark matter (CDM) models, with a cosmo- 
logical constant and a broken-scale-invariant (BSI) steplike primordial spectrum of 
adiabatic perturbations, previously found to be in very good agreement with observa- 
tions. Performing a Fisher matrix analysis, we show that in case of a large gravitational 
waves (GW) contribution some free parameters (defining the step) of our BSI model 
can be extracted with remarkable accuracy by the Planck satellite, thanks to the 
polarisation anisotropy measurements. Further, cosmological parameters can still be 
found with very good precision, despite a larger number of free parameters than in 
the simplest inflationary models. 

Key words: cosmology:theory - early Universe - cosmic microwave background. 



1 INTRODUCTION 

Current observations of the large-scale structure in the 
Universe, and of cosmic microwave background (CMB) 
anisotropics in particular, already allow a discrimination 
among different cosmological scenarios with increasing pre- 
cision. Nevertheless, many possibilities are still viable, with 
different assumptions concerning e.g. the matter content of 
the Universe, and the primordial (initial) fluctuation power 
spectrum. However, most scenarios should be excluded by 
cross-correlating the forthcoming experiments, like, for in- 
stance, balloon and satellite measurements of small scale 
CMB anisotropies, and new reshift surveys (Wang, Spergel 
& Strauss 1999). The most precise scheduled experiment 
for the measurement of the CMB anisotropies is the Planck 
satellite^], the data from which will very likely favour a re- 
stricted family of cosmological scenarios, hopefully with a 
small number of free parameters. 

As the simplest CDM model with a flat primordial 
power spectrum is already excluded, it is necessary to in- 
troduce some refinements either in the content of the Uni- 
verse (i.e., in the transfer functions of matter and radiation), 
or in the generation of initial fluctuations (i.e., in the case 



* For the instrumental specifications of the mis sion, see 
http: / /astro. estec.esa.nl/SA-general/Projects/Planck/ 



of inflationary models, in the primordial power spectrum). 
By now, the simplest variant favoured by experimental data 
seems to be that of a flat universe with a cosmological con- 
stant, fl m + Qa = 1, and a scale-invariant primordial (or 
slightly tilted) power spectrum (Kofman, Gnedin & Bah- 
call 1993; Bagla, Padmanabhan & Narlikar 1996; Ostriker 
& Steinhardt 1995; Lineweaver 1998). In two recent papers 
(Lesgourgues, Polarski & Starobinsky 1998a, 1998b, further 
referred as LPS1, LPS2), the combination of this ACDM 
scenario with an inflationary model introduced by Starobin- 
sky (1992), predicting a broken-scale-invariant (BSI) step- 
like primordial power spectrum, was investigated. In LPS1, 
the case of adiabatic primordial fluctuations was consid- 
ered when the contribution of the tensorial fluctuations to 
the CMB anisotropies is negligible. In LPS2, the possible 
contribution of gravitational waves to the CMB anisotropy 
was taken into account: it is a most interesting peculiarity 
of these models that these distinct cases are possible and 
were shown to be viable regarding observations. Using po- 
larization measurements on the precision level scheduled for 
Planck, this large GW contribution will allow accurate pa- 
rameter extraction. Briefly, the motivations for considering 
steplike models are the following. First, an even better agree- 
ment with the data can be found than in the case of a flat 
or tilted spectrum, inside a wider region of the cosmological 
parameter space. Second, a few authors point out the pos- 



© 0000 RAS 



2 J. Lesgourgues, S. Prunet and D. Polarski 



sible observational evidence for a spike in the matter power 
spectrum at fc ~ 0.05 h Mpc -1 (Einasto et al. 1997a, 1997b, 
1997c; Retzlaff et al. 1997; Gaztaiiaga & Baugh 1998). This 
is still a point of debate; however, as our BSI model predicts 
a similar feature, it would clearly be an excellent remain- 
ing candidate if the spike were to be confirmed by future 
redshift surveys. 

One could argue that BSI inflationary models, by intro- 
ducing additional free parameters in the primordial power 
spectrum, just increase the degeneracy among different sce- 
narios; that instead of making real predictions, like the sim- 
plest slow-roll models, they are just introduced ad hoc, in 
order to fit any observations; and finally, that in the case 
of BSI primordial spectra, the prospect of extracting the 
cosmological parameters at the per cent level with Planck 
would collapse. However, we recall that our model is based 
on simple assumptions concerning the inflaton potential, and 
cannot be tailored at will in order to fit any given obser- 
vational data. On the contrary, it predicts a very peculiar 
observable feature in the matter power spectrum at inter- 
mediate scales (~ 125 /i" 1 Mpc) while it makes, of course, 
predictions on all scales both for the matter and the radia- 
tion power spectrum. Even when only the radiation power 
spectrum is considered, we show in this work the following 
points: 

• the future Planck results should easily discriminate be- 
tween our BSI model and other scale-free models. 

• assuming that this model is indeed realized in nature, 
in spite of four additional degrees of freedom in the theory 
compared to the simplest versions of inflation, Planck should 
still be able to measure accurately both the cosmological and 
the inflationary parameters. Furthermore, it turns out that 
one of the inflationary parameters, p, which defines both the 
height and the shape of the step, should be constrained with 
remarkably high precision, a fact which could be of signif- 
icant interest for building particle-physics-motivated infla- 
tionary models. 



2 THE MODEL 

We assume for simplicity that our Universe is known to 
be spatially flat, and that neutrino mass and reionisation 
can both be neglected (relaxing these assumptions would of 
course increase the uncertainties on all parameters). Then, 
our model contains three cosmological parameters (h, Qa, 
fib), and five inflationary parameters, which can be under- 
stood as follows: 

• the power spectrum of adiabatic perturbations has a 
scale-invariant tilt n s on large scales, k < ko, undergoes a 
break (the shape of which is defined by one single parameter 
p) at k > ko, and is finally flat on small scale k S> ko. The 
ratio between the power spectrum on the small-scale plateau 
and at fco is given by p~ 2 . 

• the spectrum of GW has no break at fco, while the ten- 
sor tilt on small scale fc > fco is irrelevant for our purpose, 
because the corresponding contribution to the Ci's is neg- 
ligible. Using the slow-roll conditions valid on large scales, 
the scale-dependent tilt riT(k) for fc < fc can be found as a 
function of ns and nr(fco). 



• as the slow-roll approximation is still valid for large- 
scale perturbations, at fc = fco one can relate the ampli- 
tude of the GW power spectrum to the dimensionless pa- 
rameter H% Q G, and the scalar power spectrum amplitude to 
H 2 ka G/n T {k ). 

In summary, the five free inflationary parameters are: 

1. H\ G, the overall dimensionless normalization factor. 
Varying H^G (all other parameters being fixed) is ex- 
actly similar to varying the commonly used Qio, the 10- 
th multipole of the temperature anisotropy power spectrum 
(Lineweaver & Barbosa 1998). Hence, we will further use 
this parameter instead of H% G. 

2. fco, the scale of the break. 

3. p, which defines the break's amplitude and shape. 

4. ns, the scale-invariant scalar tilt on scales fc < fco. 

5. n-r(fco), the (effective) tensor tilt at fco. 

The usual tensor-to-scalar ratio Ci /Cf does not appear 
in a natural way in this description. For fixed values of the 
cosmological parameters, there is a non-trivial dependence of 
Cio/Cfo on riT(fco) and ns- Therefore, fixing the parameters 
nr{ko) and ns fixes the ratio Cf /Cf as well. 

In previous studies (LPS1, LPS2), Q b h 2 = 0.015 was 
assumed while fco was fixed by the Einasto et al. clus- 
ter data (the spike in the matter power spectrum at fc = 
0.05 h Mpc -1 requires fc = 0.016 h Mpc -1 ). Two possibil- 
ities for the scalar tilt were investigated: 

A. n s ?a 1, which implies dnr/dhik w n%. 

B. n s w 1 + rir = constant, which implies n-r(fc) = 
nr{ko) = constant. 

Further, a double normalization was performed to both 
Qio = 18 ij,K (Bennett et al. 1996) and cr 8 = 0.60 fi -0 ' 56 
(White, Efstathiou & Frenk 1993). With these constraints, 
the remaining free parameter space was three-dimensional: 
(h, Q,a, nr(fco)), or equivalently (h, SIa,CT /C¥ ). In 
both cases A and B, the preferred regions following from 
the current observations were found. Now, we choose 
one point inside each allowed region, corresponding to 
(h,n A ,C[ /Cf ) = (0.7,0.7,0.8) for both cases A and B. 
Assuming that each of these two points describes the "true" 
cosmological scenario, we perform a Fisher matrix analy- 
sis with eight free parameters (h, Qa, fib, Qio, fco, P, ns, 
n T (k )). 



3 THE FISHER MATRIX 

Using the CMB Boltzmann code CMBFAST (Seljak & 
Zaldarriaga 1996), we compute the derivative of the Ci's 
with respect to each parameter 9i t i=i,.,8. The Fisher matrix 
(Jungman et al. 1996a, 1996b; Tegmark, Taylor & Heavens 
1997; see also Bond, Efstathiou & Tegmark 1997; Copeland, 
Grivell & Liddle 1998; Stompor & Eftathiou 1998; Eisen- 
stein, Hu & Tegmark 1998) is then obtained by adding the 
derivatives, weighted by the inverse of the covariance ma- 
trix of the estimators of the polarized and unpolarized CMB 
power spectra for the Planck satellite mission , Cov(C*, Cj): 

^ = E£^ov- ( cM)^, (i) 

(=2 X,Y 
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where {X, Y} G {T,E,TE} (Kamionkowski, Kosowsky, & 
Stebbins 1997, Zaldarriaga, Spergel & Seljak 1997; Prunet, 
Sethi & Bouchet, 1998a, 1998b). The meaning of Fij is the 
following. Assuming that a fit to the Planck data yields a 
maximum likelihood for the model under consideration (for 
which the derivatives were computed), the 1-a confidence 
region in the eight-dimensional parameter space would be 
inside the ellipsoid (Press et al. 1989): 

^MiAOjFij =9.3 . (2) 

Using Fij (or the dimensionless Fisher matrix Fj = 
OiOjFij), one can also compute the allowed region in lower 
dimensional cuts of the parameter space, making no assump- 
tions on other parameters. In particular, the 1-a uncertainty 
on a single parameter is just the square root of a diagonal 
coefficient of the inverse Fisher matrix: 

A6i = ^{F-^)u . (3) 

Each multipole will be measured by Planck with a preci- 
sion of the order of 1%. As there are many more independent 
measurements than free parameters, one naively expects the 
parameter extraction to be much more precise. However, in 
general, the parameters are degenerate, i.e., some combi- 
nations of parameters produce a very weak change in the 
Ci curve. Hence, even when some other combinations can 
be measured with very high precision, each parameter sep- 
arately is constrained only at the percent level (unless its 
effect is "orthogonal" to the other ones). 

A usefull way to express the results of the Fisher ma- 
trix analysis, which does not depend on a particular choice 
of basis in the parameter space, and contains the most re- 
fined constraints that can be deduced from the experiment, 
is to diagonalize Fij. The eigenvectors correspond to the 
axes of the likelihood ellipsoid, and the inverse square root 
of the eigenvalues to the 1-a relative uncertainties on each 
eigenvector. Eigenvectors with the smallest uncertainties are 
the best constrained parameter combinations (they generate 
maximal changes in the anisotropy curve) . Eigenvectors with 
the largest uncertainties are the worst-constrained combi- 
nations (they generate minimal changes), and are generally 
called degenerate directions in parameter space. 

In computing the covariance matrix of the CMB power 
spectra, we accounted for the presence of foregrounds (both 
polarized and unpolarized) in the measurement of the CMB 
power spectra, using the method described in Bouchet, 
Prunet & Sethi 1998 (see also Prunet, Sethi & Bouchet 
1998a, 1998b). The derivatives of the spectra have been com- 
puted with an error-minimizing routine derived from Nu- 
merical Recipies (Press et al. 1989). 



4 RESULTS 

The uncertainty on each parameter is presented in Table 1. 
A and B stand for the two models previously mentioned. 
T is a tilted model, with only three inflationary parameters 
instead of five, namely (Qio, ns, nr(fco)), with the same 
values as in model B. Of course, one should keep in mind 
that all the uncertainties quoted in Table 1 would increase if 
the space of free cosmological parameters were enlarged. In 



the lines without a x-sign and with italic numbers, the po- 
larisation measurement is not taken into account. The main 
conclusions to be drawn from the table are the following. 

First, the three cosmological parameters are constrained 
with almost the same precision for the tilted and for the BSI 
models; this means that the step parameters (fco, p) do not 
"conspire" with the parameters (h, Qa, tlb) hi order to create 
directions of degeneracy. Hence, in general, the one percent 
parameter extraction proposed by Planck is not affected in 
the case of BSI steplike models. 

The situation is somewhat different for the inflation- 
ary parameters. The normalization and tilts, (Qio, ns, nr), 
appear less constrained; on the other hand, the step param- 
eters, (ko, p), can be predicted with excellent accuracy, up 
to a 0.09 % 1-a errorbar for p ! These results can be easily 
understood, especially if one keeps in mind that the best 
constraints come from high I multipoles, for which the cos- 
mic variance can be neglected. For tilted models, the scalar 
tilt enters in all multipoles, and can be accurately deter- 
mined from high l's; the two remaining inflationary param- 
eters (Qio, nr) have a similar effect on high l's (since nr 
is proportional to the tensor to scalar ratio), but measure- 
ments of the Ci's for small Ps and polarisation measurements 
reduce the degeneracy. For our BSI models, the scalar tilt 
cannot be deduced from high l's (it is defined at k < ko, i.e., 
mainly I < 100); the three parameters (Qio, ns, nr{ko)) 
combine into several degeneracy directions that can be re- 
solved only by small I measurements, so the precision re- 
mains poor. The situation is exactly opposite for the step 
parameters (fco, p), which have the crucial property of play- 
ing a role only at I > 150. Hence, they are only marginally 
affected by cosmic variance. Further, p is orthogonal to the 
degeneracy directions, and can be extracted with great pre- 
cision. 

All these features can be deduced with more accuracy 
from the Fisher matrix diagonalization, given in Table 2. 
The first lines give parameter combinations that are con- 
strained with great precision; the last lines indicate the di- 
rections of degeneracy in parameter space. It is straightfor- 
ward to see that the inflationary parameter p contributes 
only to the first four lines. Therefore, it doesn't suffer from 
any degeneracy, and is the best constrained parameter. The 
5th and 6th lines show that a change in fco (resulting in 
a slight change in the location of the first acoustic peak, 
through the change in the primordial power spectrum), can 
be cancelled by a change in the cosmological parameters 
(i.e., in the sound horizon scale). The 6th eigenvector has a 
1.5% uncertainty: this is already a small degeneracy, and h, 
Qa, £lb, ko are not as well constrained as p. The last two lines 
show the large degeneracy between Qio, ns and nr(/>o). 

Let us compare the uncertainties when polarisation 
measurements are taken into account and when they are 
not. Usually, the addition of polarized spectra leads to a 
small precision increase (by a factor < 2) for all parame- 
ters which were not part of a specific degeneracy, as can be 
seen e.g. for tilted models, by comparing the (T) and (Tx) 
data in Table 1. In our BSI model, the precision on (ns, 
nr(ko), p) increases by a much larger factor, and even by 
one order of magnitude for the parameter p ! So, measuring 
the polarisation is even more important when one consid- 
ers primordial spectra with additional free parameters (i.e., 
additional potential degeneracies to remove). One could be 
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Table 1. In the upper part of the table, we give the parameter values for the chosen mod- 
els: two BSI models A and B with 8 free parameters, and one tilted model T, with 6 free 
parameters. We also indicate the related value of C[ /Cf . The corresponding relative 1-cr 
uncertainties, A0,/#j, are given in the lower part, in percent. In the lines without a x-sign 
and with italic numbers, the polarisation measurement is not taken into account. The uncer- 
tainty on C^j/Cj'q was not calculated, but it is of the same order as the one on ny(fco) since 
in a first order description, C[ /Cf is approximately proportional to jvr'(fco). 
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Table 2. Orthonormal eigenvectors of the dimensionless Fisher matrix Fij, with their 1- 
a uncertainty (in percent). The first lines show some combinations of the parameters that 
can be recovered with a precision much smaller than 1 %. The last lines correspond to the 
directions of degeneracy in parameter space. 
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surprised by the factor 10 found for p in model B. In fact, 
when polarisation is not taken into account, p enters into a 
single combination of parameters leading to a degeneracy. 
When polarisation is added, this degeneracy is supressed 
and, as we saw, p doesn't enter into any degeneracy at all. 
This mechanism is illustrated in figure 1. 



5 CONCLUSION 

In this paper, we considered an inflationary model with BSI 
primordial spectrum and we investigated the precision with 
which the cosmological parameters and the free inflation- 



ary parameters could be extracted by the Planck satellite. 
We first conclude that in the framework of the BSI steplike 
models considered here, the extraction of cosmological pa- 
rameters can be as precise as in the case of tilted models. 
The step parameters p and fco can be constrained with ex- 
cellent accuracy, especially p, the effect of which on the C\ 's 
can be easily distinguished from the effect of any parameter 
combinations. There is no degeneracy with tilted models, 
which are special cases of our model with respect to the 
CMB anisotropies whenever ko > 0.25 h Mpc -1 . Further, 
if this class of models (or some other BSI model) were ever 
confirmed by future observations, it would be reasonnable to 
expect constraints on some of the inflaton Lagragian param- 
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Figure 1. 1-<t likelihood regions in the (Ap/p, A!Ja/^a) plane, 
with and without including polarisation measurement. The only 
degeneracy involving p is removed by the introduction of polar- 
ization measurement. Therefore, the ellipse appears vertical in all 
(Ap/p, A9i/6i) plots. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
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eters up to the 0.1% precision level. This is most interesting 
for building particle physics inspired inflationary models. 
On the other hand, precision is lost for the determination 
of the scalar and tensor tilts on large scales, as well as on 
the quantity Cx /Cf , related by the slow-roll equations to 
nr(fc). Finally, in usual inflationary models, the inclusion of 
polarization measurements is known to increase the preci- 
sion for the parameter extraction. In our model, polarization 
measurements by Planck are even shown to render the ex- 
traction of the inflationary parameters up to about 10 times 
more accurate. 
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